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^O . A Ginsparg- Wilson based calibration of the topological charge is used to calculate 

^— V I the renormalization constants which appear in the field-theoretical determination of 

CN I the topological susceptibility on the lattice. A systematic comparison is made with 

calculations based on cooling. The two methods agree within present statistical 
Q ■ errors (3%-4%). We also discuss the independence of the multiplicative renormal- 
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ization constant Z from the background topological charge used to determine it. 
PACS numbers: ll.15.Ha, ll.lO.Gh 



INTRODUCTION 



Topology in QCD plays a relevant role in understanding several low energy properties 

of the theory. One of the quantities that have a direct interest in phenomenology is the 

c^ . topological susceptibility x which is defined as the correlator at zero momentum of two 

topological charge density operators Q{x). In particular the value of x iii the pure gauge 

theory is interpreted in terms of the mass of the singlet pseudoscalar meson [l|, |^ . 

>-Ch ■ The lattice is an excellent tool to calculate such dimensionful observables and specifi- 

>• . cally the lattice determination of y in the pure gauge theory is in good agreement with 

K> '. phenomenological expectations [s], y, |5| . 

^ I The calculation of x or any other topology-related quantity on the lattice requires a 

regularization Ql{x) of the topological charge density operator. The formal naive continuum 

limit must satisfy Q L{x)^^a'^Q{x) , (a is the lattice spacing). 

However in general Ql{x) does not meet the continuum Ward identities 6]. Consequently 
the lattice definition of the topological susceptibility xl = {Q\) /V iY is the lattice volume 
and Ql = J2x Ql{x)) need not coincide with the physical continuum expression x- The two 
quantities are related by [T], |8| 

XL = Z'a\ + M, (1) 

where Z and M are renormalization constants which for evident reasons are called multi- 
phcative and additive respectively. In order to extract x from the lattice data of xl one 
must know the values of Z and M. 
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Let us outline the origin of the two renormahzation constants. Any matrix element 
which contains n insertions of the topological charge operator can be calculated either in the 
continuum (with an adequate regularization) or on the lattice. In general the two calculations 
match only after the inclusion of a multiplicative renormahzation constant .Z' J8|, |9| for each 
insertion of Q{x). In a formal writing, Ql{x) = Za^Q{x) where Z is a finite renormahzation 
constant. In the theory with fermions the topological charge mixes with other operators 
related to the axial anomaly jll|. This mixing induces a correction to the above described 
multiplicative renormahzation and consequently also to Eq.(^. Such a correction is however 
rather small [11] and it is usually neglected consistently with the large statistical errors from 
a typical numerical simulation. 

If n > 2 the above matrix element includes further divergences which are the origin of 
the additive term M. In fact the expression for x contains the product of two topological 
charges at the same spacetime point. The operator expansion of this product at short 
distances contains mixings with operators that share the quantum numbers of x- On the 
other hand it is known that the correlation function of two topological charge operators at 
nonzero distance is negative llj,ll3|, (Q(x)Q(0)) < for x 7^ 0. This inequality also holds 
on the lattice for any definition of Ql{x) [iJ, ll5|, lla, lUf • Since x is a positive quantity, part 
of the contact terms must be included in the physical definition of x- The rest of the terms, 
if any, must be subtracted and they are M. 

A prescription is necessary to calculate M. Due to the expression x = d^lnZ(6')/d^^|6)=o 
where Z{6) is the partition function of the gauge theory with a theta term, we know that x 
vanishes within the zero topological charge sector. We then adopt the following definition 
for M: it is the value of xl in the sector of zero topological charge, M = XL\q=o where q 
is the value of the total topological charge of a configuration (as determined by cooling or 
other means, see later). 

A nonperturbative method to calculate Z and M has been developed in Refs. J18l . [U, 
The method will be described in Section 2. It has been used in several calculations of 
topological properties in QCD and other theories. Various tests and studies of efficiency 
have also been worked out in the past. In Section 3 we will present the main contribution of 
the paper: a study of systematic errors that may affect the calculation of the renormahzation 
constants and a comparison of results when they are obtained by using different methods 
for calibrating the background topological sector. Some conclusive comments are given in 
Section 4. 



II. THE NONPERTURBATIVE CALCULATION OF M AND Z 



In Refs. J18l Ll9| a technique, called "heating method", to calculate the renormahzation 
constants in Eq.(0) was put forward. The idea behind the heating method is that the UV 
fluctuations in Ql{x), which are the ultimate cause for renormalizations, are effectively 
decoupled from the background topological signal so that, starting from a classical conflgu- 
ration of flxed topological content, it is possible, by applying a few updating (heating) steps 
at the corresponding value of the lattice bare coupling constant /?, to thermalize the UV 
fluctuations without altering the background topological content. This result is favoured by 
the fact that topological modes have very large autocorrelation times, as compared to other 
non-top ological observables (this autocorrelation time is particularly long in the case of full 
QCD [20 , 21] and also in the case of a large number of colours [22| ) . 

One can thus create samples of conflgurations within a given topological sector q with 



the UV fluctuations thermalized. If g 7^ tlien tlie measurement of Ql on such a sample 
leads to the lattice value of the background topological charge q^. As described above [8|, |9|, 
this result must be renormalized to match the continuum value g, 

qL = {QL)\, = Zq, (2) 

where the subscript \q indicates that the thermalization is achieved within the sector of 
charge q. Therefore, knowing q from the classical configuration and determining q^ from the 
measurement of Ql on the sample, Z can be extracted. 

If g = then measuring Q\/V on the sample leads to the value Xl|<?=o which is precisely 
M, as indicated above. 

M can also be calculated on samples with nonvanishing topological background q. Fol- 
lowing Eq.(^ and Eq.Q the additive constant is extracted in this case by using the relation 

M = \^{{Ql)U-{{QL)Uf) (3) 
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and leads to the same results 

A sample of configurations belonging to the sector of total topological charge q is obtained 
in the following way. One starts from a classical configuration with topological charge q. 
It can be easily obtained either by using cooled configurations where the energy and the 
topological charge correspond to the presence of one instanton (if g = 1)^ or by setting all 
gauge links to unity (if g = 0). Then a few updating steps are applied and the proper operator 
(either Ql or Q\/V) is measured at each step. Moreover at each step the background 
topological charge is checked by cooling^ j26| in order to verify that the configuration still 
lies in the sector of charge g [28] . When the result of the measurement stabilizes (data 
display a plateau), while g stays fixed, we consider that the UV fluctuations are thermalized 
and {QL)\q=i and {Q\/V)\q=Q yield Z and M respectively. 

When the cooling applied to a configuration reveals that its background topological charge 
is no longer equal to g then the configuration is discarded from the sample. However it might 
also happen the following event: a possible new instanton or anti-instanton created by the 
various updating steps might evade the cooling probe because the very cooling procedure 
could destroy it (this may happen especially when the instanton spans a few lattice spacings). 
In this case we would include in the sample a configuration which actually does not belong to 
the sector of topological charge g. Such an event would obviously distort the measurement of 
any of the above topological operators. For example, it yields an overestimation of ((5|)|q=o 
because any added instanton or anti-instanton only increments the value of the square. On 
the other hand since the theory predilects the sector of vanishing topological charge, the 
updating steps during the calculation of Z tend to bring the configuration to that sector 
either by destroying the original instanton or by creating from scratch an anti-instanton. 
As a consequence the value of {QL)\q=i becomes underestimated because on average the 
expectation value of Ql in the sector g = is less than in the sector with g = 1. 

In the past we have always been aware of that potential problem and in this paper we 
present a study where the background topological sector is calibrated by another method in 



^ Alternatively one can also approximate a BPST instanton 24] on the discrete lattice [251 ■ The same 

procedures can be applied for g > 1. 
^ Different variants of the cooling procedure lead to identical results [27| . 



order to compare the results and detect any difference in the form of a systematic error. The 
new method is the counting of zero modes by using Ginsparg- Wilson based operators [29| 
and it will be described in the next Section. 



III. A STUDY OF SYSTEMATIC ERRORS 

Following the lines described in the above Section, we have calculated the values of the 
renormalization constants in Eq.(0) for the 1-smeared topological charge operator defined 
as [33 

1 ^^ 
*5^(^) = ~^9^ Yl e^i'paTr{n^^(x)np^(x)} , (4) 



where all link matrices have been substituted by 1-smeared links [3l[ (the smearing param- 
eter was c = 0.90). The corresponding renormalization constants will be called M^^^ and 
Z^^") to denote the level of smearing. In the above expression Il^,y{x) is the plaquette in the 
ft — u plane with the four corners at x, x + fl, x + 'fl + u, x + v (counter-clockwise path). 
Links pointing to negative directions mean U^^{x) = UUx — ^). The generalized completely 
antisymmetric tensor is defined by ?i234 = 1 and 'e(^-^)upa = —^^.upa- The calculation was 
performed for the Yang-Mills theory with SU(3) gauge group and Wilson action [3^. The 
lattice size was 12^ and the bare coupling /? = 6. 

The practical procedure was the following: starting from a classical configuration with 
the adequate topological background content (g = for M^^'> or g = 1 for Z'^^}Y 80 heat-bath 
(HB) steps were applied, each step consisting of three Cabibbo-Marinari \?)S\ hits, one for 
each SU(2) subgroup and Kennedy-Pendleton algorithm ^5^ for refreshing the dynamical 
variables. The operator {{QlY /V for M^^^ or Ql for Z^^^) was measured every 4 steps. 
This set of 20 measurements is called "trajectory". A test of the topological sector was 
accomplished after each measurement on a separate copy of the configuration. We accepted 
only those measurements that were thermalized within the corresponding topological sector. 
This condition means that the data must have stabilized to a plateau and that the test must 
have revealed that the configuration lay within the correct topological sector. The average 
over all accepted measurements yielded M^^"* or Z^^\ 

The test was performed by two different methods: either by a traditional cooling !26| or 
by counting fermionic zero modes (CFZM). This last method consists in calculating the net 
number of zero modes, n_|_ — n^ by enumerating the level crossin gs i n th e sp ectrum of the 
Wilson-Dirac operator Dyv — dm as the fermion mass m is varied 35|, |3a, |37 1 . This method 
was utilized in Refs. J37|, |38[ to calculate the topological charge. 



The main advantage of the CFZM method against cooling is that the first one does not 
need to modify the configuration to which it is applied, hence the topological content is never 
altered during the test. On the contrary its main disadvantage is that its implementation is 
heavily time-consuming. 

Let us discuss in more detail the CFZM method. We can stop counting crossings at any 
mass am = arustop inside the allowed interval am < 2. Thus, in general the measurement of 
the topological charge will depend on amgtop because there can be level crossings all along 
the interval of masses where the gap is closed. This makes the zero mode counting method 
to look ambiguous. However in Ref. [37| it is shown that such a dependence is rather mild 
for /3 ~ 6 as long as arristop > 1-5. In the present work we have sought crossings up to three 
different values for the stopping mass: 1.0, 1.5 and 2.0. In particular the last (and largest 
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FIG. 1: {Q'j^/V)\q=o as a function of the HB step for the four cahbration tests described in the 
text. The results for M^^^ have been placed on the left side and are represented by white symbols. 
Data corresponding to different calibration methods have been shifted with respect to each other 
in order to render the Figure clearer. 



possible) value looks particularly interesting because instantons representing crossings close 
to am = 2.0 span a size of a few lattice spacings ^2|, i.e.: they are the instantons that most 
likely could evade the cooling test. 

The search for crossings was realized by following the same procedure of Ref. [38[. An 
accelerated conjugate gradient algorithm j^ was employed to extract the lowest eigenvalues 
of the Wilson-Dirac operator. 

The topological charge after the cooling test was required to be equal to 1 or within a 
tolerance of 6. We usually chose 6 = 0.3 although tests with other values were performed 
(see later) proving that the results are very robust against variabihty of 6. 



A. Systematic errors on M^^> 

The results for M ^^^ are shown in Fig. ^ We have calculated 1380 trajectories of 80 HB 
steps and measured the operator Q'i/V every 4 steps after checking that the background 
charge was zero. For each measurement, in the Figure we show the average over all the 
results. Actually we show the data after the IQth step because measurements after too few 
HB steps are irrelevant as the configuration is surely not yet thermalized. A plateau seems 
to set in after approximately 40 steps which indicates that thermalization has been achieved. 
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FIG. 2: Counting of trajectories that still lie in the q 
steps increases for the four calibration methods. 



topological sector as the number of HB 



Then the value of M^^^ can be extracted by averaging over the data after the AOth step. 

The error bar was calculated in the following manner: each trajectory was treated as a 
single data by averaging all accepted (i.e.: all thermalized and correctly calibrated) points 
in it. Since separate trajectories are prepared by independent Monte Carlo runs, this proce- 
dure guarantees the absence of autocorrelations. Then the average and the error are easily 
obtained from this set of 1380 data. Furthermore, due to the fact that different trajectories 
can contain a variable number of accepted points, each trajectory must be weighted by a 
factor proportional to the number of accepted points in it. 

If the number of accepted points in the trajectory t is rit and if Tj is the average over all 
accepted measurements m in this trajectory. 
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Fig. n displays the results for all four methods of calibration, cooling with S = 0.3 and 
CFZM with three different values for the stopping mass. The four white symbols on the left 
side of the Figure are the corresponding results for M^^\ They become smaller as the value 
of arristop is incremented. This effect is possibly an indication of a systematic error which 
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FIG. 3: Counting of trajectories that still lie in the q = topological sector as the number of HB 
steps increases for the cooling method with three different values of 6. 



however has little effect in the calculation since all results look compatible with each other 
within errors. It must be stressed that the statistical errors, represented by the error bars 
in the Figure, amount to about 3% which is rather small. 

In Fig. 121 we show, for each measurement, the number of trajectories for which the cal- 
ibration test gave an acceptable result, g = 0. The plot is shown for all measurements, 
thermalized or not, from the 20th HB step onward. The maximum possible number is obvi- 
ously 1380 and it falls off as the amount of HB steps is increased. The decrease is steeper for 
the CFZM with larger stopping mass. In the latest step and for the CFZM with arristop = 2 
about 30% of all trajectories have varied the topological sector, while for the cooling method 
the analogous percentage is about 17%. Such a difference still allows to obtain results for 
M^^^ that are compatible within (small) errors, as shown by the white symbols in Fig. ^ 

The values of M^^' obtained by using cooling with varying 5 are indistinguishable from 
each other (white circle in Fig. ^ for 6 ranging from 0.1 to 0.5. In Fig. El three tolerance 
parameters for cooling are compared indicating that the three tests are almost completely 
equivalent. Again the maximum possible number of trajectories with the right background 
topological charge is 1380 and it diminishes as the amount of HB steps is increased. This 
number became 1136 for 6 = 0.1 at the latest step. 
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FIG. 4: Averages of {{QL)\q) /q as a function of the HB step for the four calibration methods. 
An instanton with charge q = 0.997 was used. The values of Z^^' correspond to the height of 
the plateau and they are indicated on the left hand side of the plot with white symbols. Data 
corresponding to different calibration methods have been slightly shifted for clarity. 



B. Systematic errors on Z^^' 

Fig. 0] displays the analogous study of Fig. ^for the calculation of Z^^\ An instanton 
with topological charge (as measured with Ql after cooling) q = 0.997 was employed. The 
procedure for the calculation of Z^^^ resembles very much that of the M^^\ We prepared 840 
independent trajectories and again the calibration was performed by four different tests, as 
indicated in Fig. HI The average per trajectory is Tj, 



T, 



^T.(Q^ 



and the result for Z^^^ is 
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0.997 Zt^t 

A plateau sets in at about the 20th HB step. Figures similar to Fig. El and 01 are obtained 
with analogous conclusions. Again the coincident results (white symbols in Fig. |3|) indicate 
that all calibration methods are equivalent in such a way that any systematic error implicit 
in our method has negligible consequences within our precision (statistical errors in Fig. 01 
are about 4%). 



TABLE I: Z^^> for the l-smeared operator as a function of (3 for several values of the background 
topological charge q. Instantons with charge q « —1, ^ +1 and w +2 are used on a lattice volume 
16l 



/3 


q = +0.990 


q = -0.988 


q = +0.955 


q = -0.902 


g = +1.941 


6.00 


0.373(20) 


0.383(15) 


0.370(20) 


0.365(20) 


0.390(12) 


6.20 


0.432(6) 


- 


- 


- 


0.441(5) 


6.50 


0.503(11) 


— 


— 


- 


0.506(5) 



TABLE IL Z'^' for the 2~snieared operator as a function of /3 for several values of the background 
topological charge q. Instantons with charge g Ri — 1, « +1 and ps +2 are used on a lattice volume 
16l 



p 


q = +0.990 


q = -0.988 


q = +0.955 


q = -0.902 


g = +1.941 


6.00 


0.500(12) 


0.503(13) 


0.495(15) 


0.487(20) 


0.499(16) 


6.20 


0.560(6) 


- 


- 


- 


0.569(4) 


6.50 


0.620(10) 


- 


- 


- 


0.631(3) 



It is well-known that the total topological charge of isolated classical instantons in general 
is not equal to integer numbers when it is calculated with the operator of Eq.(0]). The 
difference between the value of the charge and the closest corresponding integer (|1 — g| in 
our case) becomes negligible when the inequalities a ^ p ^ La hold (p being the instanton 
size and L the lattice size). A simple calculation shows that the value of q for a discretized 
instanton in a volume L^ is given by 



g ^ 1 



.La) 



(9) 



while the discretization error is totally negligible if a < p/3. The ratio a/ p can be estimated 
by looking at the action density distribution. We expect that the infrared effect described 
in Eq. (jU)) cancels out if we divide the charge after heating by the initial value q as indicated 
in Eq.(j2I). This fact was carefully checked in J40| for the 2D 0(3) nonlinear sigma model. 
It was also used in Eq.® although in that case and within our errors, the q = 0.997 in the 
denominator is indistinguishable from 1. In the present study we have verified it for our 
gauge theory: in Table HI and ITTl the multiplicative constant is calculated at several values of 
the gauge bare coupling /? on a 16^ volume starting from various initial instantons for the 1- 
and 2-smeared operators (they are constructed as in Eq.(JD) after substituting all links with 
1- and 2-smeared links respectively J3l|)- A number of trajectories ranging from 200 to 500 
were used. Notice that within errors the values for Z^^^ and Z^'^^ display a dependence on (3 
but not on q as long as data are divided by the initial value q, as described above following 
Eq.Q. If instead data were divided by the integer closest to q then the results for Z^''\ 
{i = 1, 2) would display a fake dependence also on q. This fact is seen in Fig. El where the 
values of Z^"^^ are displayed as a function of q after dividing by q (squares) or by the closest 
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FIG. 5: Z^'^'{(3 = 6.0) as obtained by using four several initial classical instantons (with \q\ =0.902, 
0.955, 0.988 and 0.0990). Squares stand for the results obtained by applying Eq.Q and triangles 
by applying the same equation after rounding the value of q to its closest integer. 



integer to q (triangles): only the squares show constancy with respect to q. Within our 
statistical errors the systematic error is negligible down to g ~ 0.9. Possibly if the statistics 
were increased this limit could grow. In any case to calculate Z it is better to make use of 
instantons with q as close as possible to an integer. 

IV. CONCLUSIONS 



We have studied two possible sources of systematic errors in the nonperturbative deter- 
mination of the renormalization constants for the evaluation of the topological susceptibility 
on the lattice, Z and M (see Eq.([T))): 

• i) The cooling test of the background topological charge along the heating process 
could yield a wrong information since the cooling procedure could delete the unwanted 
instanton or anti-instanton created by the heating steps. An independent check of the 
cooling test has been performed by applying calibration methods based on the counting 
of fermionic zero modes. No sizeable systematic effects have been observed within our 
(rather small) statistical errors (3% for M and 4% for Z). 



ii) In the calculation of Z for the operator Ql of Eq. (j3|) the simulation must be started 
from a configuration with a topological charge q different from zero. In the infinite 



11 

volume limit q takes on integer values, however on the lattice, the determination of 
q by using the operator in Eq.(JH), leads to results which in general are close to but 
not strictly equal to integers. We have argued that this is a potential source of error 
that can affect the calculation of Z. However it can be kept under control if one uses 
Eq.((2I) to extract Z with q not rounded to its closest integer. 

In conclusion one can safely use the field-theoretical method to study topology on the 
lattice since any possible systematic error of the method is well under control. Moreover the 
method is much less demanding in computer time than the Ginsparg- Wilson based method. 

The APEmille facility in Pisa was used for part of the runs. 
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